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Abstract 

By using results of highly accurate computations of the total energies of a large number of few- 
electron atoms we construct a few interpolation formulas which can be used to approximate the 
total energies of bound atomic states. In our procedure the total energies of atomic states E are 
represented as a function of the electric charge of atomic nucleus Q and the total number of bound 
electrons Ng. Some general properties of the E{Q,Ne) function are investigated. The knowledge 
of the E{Q,Ne) function allows one to determine the total (and binding) energies of these states 
in arbitrary atoms and ions with different Q and N^- 
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In this short communication we discuss accurate and relatively simple interpolation for¬ 
mulas which can be used to predict the total energies and other bound state properties in 
various atoms and ions. Recently, a very substantial progress has been achieved in highly 
accurate computations of the bound states in few-electron atoms and ions. Based on the 
results of such calculations we can construct a number of different interpolation formulas 
for the total and binding energies of atomic few-electron systems, i.e. for atomic systems 
with different nuclear charges Q and different number of bounded electrons Ng. Analo¬ 
gous formulas can be constructed for other bound state properties, including inter-particle 
delta-functions, various single-, two- and three-particles properties and properties which are 
determined by the expectation values of some singular operators. In this study we shall 
not discuss interpolation formulas for arbitrary bound state properties. Furthermore, below 
we restrict ourselves to the consideration of the total energies of the bound states in the 
non-relativistic atoms and ions. First, note that the overall accuracy of the Q~^ expansions 
constructed for the ground states in the two-, three- and four-electron atomic systems can 
be considered as outstanding. In this study we want to make the following step and develop 
(and later apply) the ‘universal formula’ proposed in [l|, which allows one to predict the 
total energies of arbitrary atomic systems with the known values of Q and Ng to very high 
numerical accuracy and without actual atomic calculations. In other words, we need to 
construct the universal function E{Q, Ng) of the two (integer) arguments and investigate its 
properties. 

At the hrst step of our analysis we need to determine the total energies E and non- 
relativistic wave functions as the solutions of the Schrodinger equation |2| for the bound 
states = E^, where E < 0 and H is the non-relativistic Hamiltonian of the A^e~6lectron 
ion/atom 

*2 ATe Afe 1 Ne-l Ne 1 

H = -r^TVj-Qe^T — + e^ y: V - (1) 


2m, 


e i=l 


i=l An 


i=l i=2{i>j) o 

where Vj = f = 1, 2,..., Ag and the notation n stands for the atomic nucleus. 

In Eq.([T]) the notation h designates for the reduced Planck constant, i.e. h = ^, and e is 
the elementary electric charge. In this study all masses of the atomic nuclei are assumed 
to be inhnite. In general, it is very convenient to perform all bound state calculations in 
atomic units where h = l,mg = 1 and e = 1. In these units the velocity of light in vacuum 
c numerically coincides with the inverse value of the dimensionless hne structure constant. 
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i.e. c = a where a = ~ 7.2973525698-10 In atomic units the Hamiltonian, Eq.([T]), 

is written in the form 


1 Ne. Ne 1 Ne — 1 Ne 1 

e ^ (2) 

^ i=l i=l i=l i=2{i>j) 

As follows from Eq.(l2]) the Hamiltonian 77 is a continuous operator-function of the nu¬ 
clear charge Q, or in other words, the nuclear charge Q is the continuous parameter (or 
control parameter) of this Hamiltonian. By applying the Poincare theorem one hnds that 
all eigenvalues of the Hamiltonian 77, Eq.(j2]), and all its eigenfunctions are the continuous 
functions of Q. In applications to real atoms the nuclear charge Q expressed in atomic units 
is always an integer number. Let us assume that we know the actual atomic wave function 
d' (or I d')). Then we can reduce Eq.([T]) by the following equation for the three expectation 
values 

E={H) = -ijv,{v;) - N,Q{—) + ) (3) 

^ ^In ^ ^12 

where E is the total energy of the bound state, while (V^), (^) = (^) = are 

the expectation values of the electron kinetic energy, electron-nucleus (attractive) potential 
energy and electron-electron repulsion, respectively. In derivation of Eq. ([3]) the fact that all 
atomic electrons identical particles. As it can be seen from Eq.(|3]) the total energy 77 is a 
function of the two parameters Q and which are both integer. Formally, to determine the 
exact (or analytical) form of the 77(Q, Nf.) function we need to determine three expectation 
values mentioned in Eq.([3]), i.e. the (V^) = — (Pi) , {^In) and W)- In reality, for Coulomb 
systems one also hnds an additional condition which is widely known as the ‘virial theorem’. 
The virial is written in the form 2(T) = —(H), where T is the operator of kinetic energy, 
while V is the operator of the potential energy. The explicit forms of these operators are: 

1 A^e A^e 1 A^e ~ 1 A^e 1 

and V = -QE^+E E ^ ( 4 ) 

^ i=i i=i j=i i=2{i>j) 


Note that the virial theorem can be written in one of the following forms: E = —{T), |(H) = 
77, etc. In general, by applying the virial theorem we can reduce the total number of 
‘unknown’ expectation values in the right-hand side of Eq.(|3]) from three to two. For instance, 
the expression of the total energy 77 in terms of (^) and (^) expectation values is 


77 = (77) 


^ rin 4 ri2 


( 5 ) 
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As follows from Eq.Q to obtain the explicit formula for the E{Q, Ne) function we need 
to derive analogous formulas for the (—) and (—) expectation values. Briefly, this means 
that these expectation values must be expressed as explicit functions of Q and Nf,. In 
reality, it is not possible to derive any closed analytical expression for the (^) and/or 
(^) expectation values written as a function of Q and Ng. An obvious exclusion is the 
Thomas-Fermi method 3j. By using other methods which are more accurate than Thomas- 
Fermi method it is impossible to derive the closed analytical formulas for the (^) and (^) 
expectation values. This means that it is impossible to obtain any closed analytical formula 
for the E{Q,Ng) function. However, the ‘atomic function’ E{Q,Nf,) can be approximated 
to very good accuracy by using results of highly accurate numerical calculations for a large 
number of a few-electron atoms/ions. Formally, if the total number of bounded electrons 
Ng is hxed, then we are dealing with the so-called Q~^ expansions for the total energies of a 
number of atoms/ions with different Q. For two-electron atoms and ions such series are well 
known since the middle of 1930’s and hrst papers by Hylleraas for two-electron ions (see, 
e.g., B-0 and references therein). Analogous series for three- and four-electron atomic 
systems were not used in applications, since they provided a very modest overall accuracy 
which was not sufficient for accurate evaluations. 

At this moment the situation with accurate numerical calculations of the three- and 
four-electron atoms and ions has changed. Currently, we have a large number of highly 
accurate results for the two-electron ions and quite a few different sets of accurate numerical 
...... oMa..ed f„. th.e. and fon.-e.eCon at„n.s and .ons (see. e.g.. Q. fl and .efe^ences 

therein). These results allow us to construct some accurate Q~^ expansions which also 
describe the total energies of the three- and four-electron atomic systems. Based on these 
formulas for the Q~^ expansions for two-, three- and four-electron atoms and ions we can 
try to guess the explicit formulas for an universal function E{Q, Ne), where Ng is the total 
number of bounded electrons. This problem has a fundamental value for whole atomic 
physics. Indeed, if we know the E{Q, A/) function, then we can predict the total energy of 
an arbitrary atom/ion with the given Q and Ng to high accuracy which is sufficient for many 
actual problems known from stellar astrophysics, physics of high-temperature plasmas, etc. 
In reality, accurate predictions of the total and binding energies of atoms/ions with different 
Q and Ng was an ultimate goal for many generations of atomic physicists. In this study we 
show that currently we are very close to fulhll this goal. 
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Accurate non-relativistic energies i? of a large number of the ground states in the two-, 
three- and four-electron atoms/ions can be found in Table I. All total energies are given in 
Table I in atomic units. The total energies of the two-electron (or helium-like) atoms and ions 
from Table I have been determined to very high accuracy with the use of our computational 
procedure which allows one to perform calculations with 3500 exponential basis functions 
in the wave functions. All basis functions are written in the relative coordinates r^ 2 ,f' 3 i and 
r 2 i. For three-electron (or lithium-like) ions (and Li atom) such energies have been taken 
from 4], while for four-electron atoms/ions they were chosen from j^. The overall accuracy 
of the ground state energies in three- and four-electron atomic systems is still signihcantly 
lower than the analogous accuracy for two-electron atoms and ions (see Table I). Note also 
that the Li“ ion, i.e. atomic system with Q = 3 and = 4, is bound, i.e. its ground 
2^S'-state is stable, but our current variational results for this system {e.g.,E = -7.5007185 
a.u.) are not highly accurate. This is the reason why we exclude this ion from Table I. By 
using the total energies of all atoms/ions mentioned in Table I (with the same value of N^) 
one can determine some numerical coefficients in the Q~^ expansion 


E{Q) ~ 0 ‘ 2 Q^ + ciiQ -l- 0-0 + biQ ^ biQ ^ -|- 62Q ^ ^ -f-... 


( 6 ) 


where 02 , 01,09 are the coefficients of the regular part of the Laurent expansion (or series), 
while &i, 62 , • • • are the coefficients of the principal part of the Laurent series E{Q), Eq.(l 6 ]). 
Note that the Q“^-expansion (or Q“^-series), Eq.([ 6 ]), is a typical ‘asymptotic expansion’. 
Briefly, this means that after some n > Umax all coefficients in Eq.([ 6 ]) rapidly increase 
with n. Contributions of the corresponding terms also rapidly increase with n and the total 
sum computed with the use of Eq.(j 6 ]), which includes such ‘growing terms’, has nothing 
to do with the original problem. Briefly, this means that in such cases the Q~^ expansion 
cannot be used to approximate the actual total energies. To avoid this problem we need to 
restrict the total number of terms in Eq.([ 6 ]). For instance, if we use = 18 values of the 
total energies E{Q) (computed for eighteen different values of Q), then the total number of 
terms in Eq.(| 6 ]) can be 8 , or 10, but not 16, or 18. An universal criterion can be formulated 
in the following form: overall contribution to the total sum, Eq.(l 6 ]), from the last term must 
be smaller (and even much smaller) than analogous contribution from the pre-last term. 

First ten coefficients of the Q~^ expansion, Eq.([n]), determined from the results of highly 
accurate numerical calculations of few-electron atoms and ions mentioned in Table I can 
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be found in Table II. These coefficients can be used in applications of the Q~^ expansion, 
Eq.(l6]), to other two-, three- and four-electron atoms and ions. The overall accuracy of the 
Q~^ expansion, Eq.([n]), for the total energies of the ground states in these atoms and ions is 
outstanding and can be evaluated as 1-10“®—5T0“^^ a.u. of the total energies. In general, the 
Q~^ expansion, Eq.([S]), can be used in applications to different iso-electron atomic systems. 
However, due to numerous problems in accurate computations of the three-, four- and many- 
electron atoms and ions the most successful applications of the Q~^ expansion are still 
restricted to the two-electron (or helium-like) atoms and ions. It should be mentioned 
that similar Q~^ expansions can be used for other bound state properties, e.g., to predict 
interparticle distances, expectation values of some delta-functions, etc. As mentioned above 
in this study we restrict ourselves to the total energies only and (brief discussion of the 


interpolation formulas for other bound state properties can be found, e.g., in |l|). 

Based on the results from Tables I and II we can make the new step which must lead 
to better understanding of the structure of bound state spectra in atoms and ions. Instead 
of dealing with many different E{Q) functions constructed for each series of iso-electron 
atomic systems, i.e. atoms/ions with the same Ne, we introduce a ‘universal function’ 
E{Q, Nf,) which depends upon two integer numbers: Q (nuclear electric charge) and 
(total number of bound electrons). If is hxed, e.g., Ng = 2, then the corresponding 
function E{Q,Ng = 2) = E{Q,2) must coincide with the function E{Q) known for the 
two-electron atomic systems. There are few possible approaches which can be used to 
determine the universal E{Q, Ng) function. In this study we apply the so-called direct 
approach which is based on the results presented in Table II. In fact, we shall assume below 
that coefficients presented in all columns of Table II correspond to one ‘universal’ function 
E{Q, Ng). Numerical differences in these coefficients can only be related with the variations 
in Ng. For instance, consider the hrst coefficients 02 from Table II. The exact value of these 
coefficients are —1, —| and —| for Ng = 2, 3, and 4, respectively. Therefore, the following 
general formula can be written in the form 

02 = 


^8 -f- Ng — 2^ 

(Q -\- Ng-^ 

^ 8 ^ 

^ 8 ^ 


(7) 

where Ng is the total number of bounded electrons in the atom/ion. As follows from Table 
II the analogous expression for the second coefficient can written in the following form 


«i = ^ + a^i\Ng-2) + aY’{Ng-2) 


( 2 ), 


( 8 ) 
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where and are the two unknown coefficients which are determined with the use of 
numerical values for this coefficient from the second and third columns of Table II. The 
formula, Eq.dH]), can be re-written into a slightly different form 




af - 


8 


{N, - 2) + 




(9) 


where a^i'^ and o!i^ are the coefficients from Table I for the three- and four-electron atomic 
system, respectively. The unknown value of the Oi coefficient for atomic system with 
A^e—electrons is designated in the left-hand side of Eq.(l9]) as The formula, Eq.(l9]), 

can be generalized to more complex cases, e.g.. 


AN.) 


= a, 


( 2 ) 


+ 


,(3) .(2) 


2-0 


2-0 


{Ne - 2 ) + 


.(4) 


2an ^ -|- (2n 


(2)1 (Ne-2)(N,-3) 


( 10 ) 


where is the Og coefficient from Eq.(|6]) defined for atomic systems with Nf, bound 

electrons. Analogously, for coefficients from Eq.dH]) one hnds 


+ 


l{3) .(2) 


(iVe - 2) + 


- 26l^) + 


n(iVe-2)(iVe-3) 


( 11 ) 


where notation stands for the bk coefficients {k = 1,2,3,...) from Eq.dHj) defined for 
atomic systems with bound electrons, where N^. = 2,3,4,.... This expression is, in 
fact, the Taylor-Maclaurin expansion for the coefficients upon the total number of 

bound electrons N^. Note that the formulas Eqs.(IZ|) - (fTTj) for the coefficient a 2 — ag and bk 
(A; = 1,2,...) are essentially exact upon for all two-, three- and four-atoms/ions presented 
in Table I. This means that the total energies of these ions are reproduced with the same 
accuracy which is provided by their Q~^ expansion. 

As follows from Eq.(|7|) the coefficient 02 is the ratio of two small integer numbers, and 
numerator of this fraction is a linear function of Ne. This follows from the fact that the 
main contribution to the total energy of any atom/ion comes from the which 

represents the leading term in the electron-nucleus attraction. In atomic units this term 
equals to an integer, or semi-integer number (for an arbitrary atom/ion). The coefficient 
ai for two-electron atomic systems is also a simple fraction, e.g., | for the ground states. 
The same conclusion is true for other bound states, e.g., for the 2^5'—triplet states, in the 
two-electron atoms/ions 14 1. 

In general, the analysis of the E{Q,Nf,) function(s) is very similar to operations with 
the Weizacker mass formula in nuclear physics. Such a similarity follows from the fact that 
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in both cases we are dealing with the hnite Fermi systems j^. Note that atoms and ions 
were considered as the hnite Fermi systems in atomic physics since hrst papers published 
in earlier 1960’s and [n| (see also 12| and references therein). More recent references 


can be found, e.g., in [^|. This general theory allows one to obtain the asymptotic form 
of the E{Q,Ng) function at very large Q and 12j. A few other advantages of that 
theory include the correct theoretical expressions for the ao and ai coefficients (see, Eqs.([7]) 
- dH])) which either coincide with ours (see Table II), or very close to them. The coefficients 
0 - 2 ) ^ 1 , & 2 , • • • predicted in the general theory are known only approximately, even in the case 
of the ground (bound) atomic states. For excited atomic states the accuracy of predictions 
of this theory is substantilally lower. In many cases applications of the general theory to the 
excited atomic states lead to very inaccurate predictions of ionization potentials and 
other properties. Such a situation with excited states lead to a conclusion that applications 
of the general theory to the ground states in atomic systems only. However, in this case 
we have to face the following crucial question: why do we need to perform dozens of highly 
accurate, atomic computations for construction of the accurate interpolation formula for the 
total energies, if in modern atomic physics highly accurate computations of any ground state 
in arbitrary atom/ion are significantly faster! It is clear that if we cannot generalize our 
interpolation formulas to the excited atomic states, then chances of the ‘general theory’ 101 


12| to survive in the future are very low. 


Our approach based on the use of highly accurate computational data for different 
atoms/ions allows one to reconstruct the E{Q,Ne) function can be applied, in principle 
to any atomic state, including excited states. As is well known from atomic spectroscopy, 
the bound state spectrum of any multi-electron atom/ion is represented as a set of different 
terms, where each term has its unique quantum numbers of angular momentum L and to¬ 
tal electron spin S. For our analysis this means that the function E{Q, Nf) (total energy) 
dehned above must be labeled by the two indexes L and S, which are good (= conserving) 
quantum numbers for an isolated atom/ion with bound electrons. These two quantum 
numbers are the labels of the corresponding atomic term. The bound state spectrum of 
any atom is represented as a combination of different LF-terms. An ultimate goal is to 
approximate the total (non-relativistic) energies of all bound states from all p ossible terms 


which can be found in real atoms and ions. In the ‘general theory’ [10| and 


ll| this problem 


is extremely complex, since often the atomic LS'—term does not exist in atoms/ions with 
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fewer electrons. Formally, this means that we need to construct the Q~^ expansions for each 
different atomic LS*—term. 

In part, such a strategy already works for the ground atomic states. Indeed, the ground 
states in different atoms correspond to the different terms, e.g., the ground state in the B- 
atom (hve bound electrons) is the 2^P-state, analogous state in the carbon atom is the 2^P- 
state. It is clear that by including future highly results for atoms/ions with larger number of 
electrons we can reach the ground states of the Fe atom (®D 4 -term), Co atom 9 -term), 
etc. However, if it is possible to connect all ground atomic states by one interpolation 
formula, then we can do the same for the excited atomic states too. Let us describe our 
current vision of this two-stage procedure. At the hrst stage the approach is based on the 
formulas, Eqs.(I7]), dH]), flTU]) and flTT]) allows one to determine other unknown coefficients 
ao, ai, 02 , &i, 62 ,.... The arising formulas are simple and convenient in applications to atoms 
and ions. The known coefficients Oq, Oi, 02 , 61 , 62 ,... are used in the ‘usual’ (^“^-expansion 

E{Q, N,) = a2(iVe)Q' + a,{Ne)Q + ao(iV,) + b,{N,)Q-^ + h{N,)Q-^ 

+ b2{Nf,)Q ^ -|- h^{Nf.)Q ^ . (12) 

where now all coefficients are the functions of the total number of bound electrons N^.. In 
this study all coefficients in Eq. ffT 2 l) are constructed as polynomial functions of Ne. The 
(5“^-expansion, Eq.([T 2 ]), provides high numerical accuracy for the total energies of all two-, 
three- and four-electron atomic systems presented in Table I. Very likely, that the analytical 
expression for each of the hi{Nf.) coefficients (i = 1 , 2 , ...) in Eq. flT^ is more complicated 
than a simple polynomial in N^.. It is clear that such an expansion must also include the 
negative powers of A),. In reality, we can investigate the E{Q,Ne) function in detail when 
we can obtain highly accurate results for hve- and six-electron atoms and ions. When highly 
accurate computations of the hve- and six-electron atoms will be completed, then we can 
substantially improve our current knowledge of the E{Q, Nf.) function(s). Moreover, we can 
derive some compact and accurate formulas for numerical approximations of these functions 
for diherent bound atomic states. This is an answer to an old question about possibility to 
hnd an analytical formula for the total non-relativistic energies of bound states in atomic 
system which contains Ng bound electrons moving in the held of inhnitely heavy nucleus 
with the electric charge Qe. In this study to approximate the E{Q, Nf.) function we have 
restricted to the ten-term formula for the Q~^ expansion, Eq.(l 6 ]). In part, such a restriction 
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is related with relatively low accuracy of the computational energies obtained for the three- 
and four-electron atoms/ions. 

At the second stage of the procedure we need to find relations between coefficients of 
these series, dehne these coefficients as the functions of etc. As mentioned above some 
bound LS”—state (or LS'-term) may not exist for atoms/ions with fewer electrons. Therefore, 
we need to solve the problem of genealogical relation between terms (or bound sates) in 
atoms/ions with different number(s) of bound electrons Ne. To explain this problem let us 
assume that we have determined known function Els{Q, Ne) (total energy) for atoms/ions 
with the same number bound electrons Ng {Q is varied). Now, suppose that the total number 
of bound electrons increases by one, i.e. A/ —)■ A'e + 1- Now, we have the new term L'S' 
and the new function El/s'{Q, Ne + 1) (total energy). We need to predict possible numerical 
values of L' and S' quantum numbers. As it follows from the fundamental principles of atomic 
theory for the new spin quantum number we have ^" = ^±1, ifS' 7 ^ 0 , and S' = S + ^, 
if S' = 0. To predict the new value of L' we have to know the angular momentum i of the 
additional electron. If we know this value, then one hnds that all possible L' values are 
located between the two following limits \ L — i \ (lower limit) and L + £ (upper limit), i.e. 

\ L — £ \< L' < L + £. If these conditions for the L' and S' are obeyed, then the total energies 
Els{Q, Ne) and Ei'S'iQ, Ne + 1 ) can be used in one series. 

As follows from the results of our study the function E{Q,Ne) can be constructed for 
all bound states in multi-electron atomic systems, if (and only if) their terms are directly 
connected by the ‘selection’ rule mentioned above. For instance, in this study we discuss 
the total energies of two-, three- and four-electron atoms/ions which have their ground 
1^5'—,2^5'— and 2^5'—states, respectively. These bound states have different multiplicities, 
i.e. they are singlets and doublets states. Suppose we want to include in our analysis hve- 
electron atoms/ions (B-like ions) for which the ground state is the 2^P—state. The ‘selection 
rules’ mentioned above works in this case. Therefore, we can construct the E{Q,Ne) func¬ 
tion for the , 2^5'—, 2^5'— and 2 ^P—states in the two-, three-, four- and hve-electron 
atoms/ions, respectively. However, if we want to evaluate, e.g., the total energies of the 
states in the hve-electron ions, then we have to use a slightly different consequence of 
bound states the two-, three- and four-electron atoms/ions, respectively. A natural choice in 
this case is to consider the bound ,2^5'— and 2 ^P—states in the two-, three- and four- 
electron atoms/ions, respectively. The total energies determined for the bound 3^D—states 
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in the five-electron atoms/ions perfectly complete the data (total energies) computed for the 
bound ,2^5'— and 2^P—states in the two-, three- and four-electron atoms/ions. On 
the other hand, it is clear that our total energies from the third column of Table I (total 
energies of the 2^5'—states in the four-electron atoms/ions) are useless in this case. In other 
words, these energies cannot be used to predict the total energies of the bound 3‘^D—states 
in five-electron atoms/ions. Formally, this means that to construct highly accurate inter¬ 
polation formula for the total energies of the bound 3^D—states in the five-electron atoms 
and ions one needs to perform a number of separate, highly accurate computations of the 
bound 2^P—states in the four-electron atomic systems. In general, the total energies of any 
bound state in multi-electron atomic system can be predicted to relatively high accuracy, 
if we know a number of total energies of bound states in atomic systems with fewer elec¬ 
trons. SUch bound states in atomic systems with fewer electrons correspond to the different 
//S'—terms. However, all these bound states (or terms) must be related to each other by the 
selection rules mentioned above. An additional problem follows form the fact that some of 
the atomic terms can be connected by using different atomic terms in the systems with fewer 
bounded electrons. It is clear that the procedure must be self-correlated. This means that 
we must have some additional relations between the E^siQ: A'e) and Eys'^Q^ A"') functions 
constructed for different atomic terms and for atoms/ions which contain different number(s) 
of bound electrons. This interesting question cannot be answered at the current level of 
theoretical and computational development. 

Appendix 


The formulas Eqs.([7]), ([9]), ffTOj) and flTT]) for the coefficients 02 , oi, oq and 61 , 62 , ■ • • from the 
main text can be re-written in a number of different forms which are often more convenient in 
applications. In this Appendix we present the explicit formulas for the coefficients 02 , Oi, Oq 
and 61 , 62 , • • • for five- and six-electron atomic systems (i.e. toms and ions). In all these cases 
the formula for the 02 coefficient coincides with Eq.(I71) (here we do not want to repeat it). 
The formulas for other coefficients from Eq. fll2l) take the following forms. For Ng = 5 one 
finds: 


(iVe = 5) 
“1 


+ 


,(3) 5.,^^ r(4) ,J3) , 5T(iV,-2)(iVe-3) 


a - 


(iVe-2)+ ar-2a\^> + 


8 J 



3aS^^ + 3af ^ 


5n(iVe-2)(iVe-3)(iV,-4) 

sJ 6 


(13) 
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a, 


(Afe = 5) 
0 


= + 


.(3) ^(2) 


^0 


^0 


(iVe - 2) + 


„(4) o„(3) I (2)1 - 2)(iVe - 3) 


02^^ — 302*^^ + 3a2^^ — a2^^ 


(7V,-2)(iVe-3)(iVe-4) 

6 


( 14 ) 


and 


,{N^=5) 


= if’ + 


(,(3) .(2) 


(iV, - 2) + 


^(4) 




(iVe-2)(iV,-3) 


- 36i^) + 3bf - 


(iVe-2)(iVe-3)(iV,-4) 

6 


(15) 


where A^e = 5 in the right-hand sides of these equations. 
Analogous formulas for the A'e = 6 are 


at 


{N^=e) 


— t: + 


,(3) 


5 


- o 


(N,-2)+ at'-2ai"' + J 


.») 


5, (iV, - 2)(]V, - 3) 


ar-=»> = af + 


af > - 3at> + 3a? - | ^-2)^-3)^-^) 

„<«) _ 4aS''> + 6a? - 4a? + 5l (iV. - 2)(iV. - 3)(A^. - 4)(iV. - 5) 

8-1 24 

(iVe-2)(iVe-3) 


(3) (2) 

ay - ay 


(iVe - 2) + 


— 2ao^^ + a^^ 


a!? - 3a!? + 3a? - a?' (N.-2m-3m-4) 




4a?) + ha?) - 4af + a® 


(16) 


6 

(iVe-2)(iVe-3)(iVe-4)(iVe-5) 

24 


(17) 


and 


,(7Ve=6) 


= b?+ 


r(3) ,(2) 


(iVe - 2) + 


t - 26f + 1 - 3) 


4® - 34i*> + 3b? - - 2me - 3m. - 4) 

b? - 44f + 64^ - 4b? + b?] W-2)W-3)W-4)(iV.-5) 


(18) 


where in the right-hand sides of these equations we must put Ng = 6. Generalization of 
these formulas to the cases when A'e > 7 is straightforward and relatively simple. 

As follows from these formulas the addition of one electron, i.e. Ae — 1 ^ Ng, affects only 
the last term in each of these formulas, while all previous terms contains coefficients known 
for atomic systems with Ae — 1 electrons. This allows one to consider different atomic 
LA-terms for the bound states with different Ng (it is assumed that the ‘selection rules’ 
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mentioned in the main text are obeyed for these LS'-terms). 
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TABLE I: The total non-relativistic energies E of the different atoms/ions in their ground states 
in atomic units. All nuclear masses are infinite. Q is the nuclear electric charge and N^, is the total 
number of bounded electrons. 


Q 

Ne = 2 

Ne = 3 

iVe =4 


1 

-0.5277510165443771965925 



2 

-2.90372437703411959831115924519440 


-- 

3 

-7.27991341266930596491875 

-7.4780603236503 


4 

-13.65556623842358670208051 

-14.3247631764654 

-14.667356407951 

5 

-22.03097158024278154165469 

-23.424605720957 

-24.348884381902 

6 

-32.40624660189853031055685 

-34.775511275626 

-36.534852285202 

7 

-44.781445148772704645183 

-48.376898319137 

-51.222712616143 

8 

-59.156595122757925558542 

-64.228542082701 

-68.411541657589 

9 

-75.531712363959491104856 

-82.330338097298 

-88.100927676354 

10 

-93.906806515037549421417 

-102.682231482398 

-110.290661070069 

11 

-114.28188377607272189582 

-125.2841907536473 

-134.980624604257 

12 

-136.65694831264692990427 

-150.1361966044594 

-162.170747906692 

13 

-161.03200302605835987252 

-177.238236559961 

-191.860986338262 

14 

-187.40704999866292631487 

-206.5903022122780 

-224.051310298012 

15 

-215.78209076353716023462 

-238.1923876941461 

-258.741699427160 

16 

-246.15712647425473932009 

-272.0444887900725 

-295.932139288646 

17 

-278.53215801540009570337 

-308.1466023952556 

-335.622619375075 

18 

-312.90718607661114879880 

-346.4987261736714 

-377.813131866050 

19 

-349.28221120345316700447 

-387.1008583345610 

-422.503670826658 

20 

-387.65723383315855621790 

-429.9529974827626 

-469.694231675265 

21 

-428.03225432023469116264 

-475.0551425155 

-519.384810821074 

22 

-470.40727295513838395930 

-522.4072925498 

-571.575405411671 

23 

-514.78228997811177388135 

-572.0094468708 

-626.266013153662 

24 

-561.15730558958127234352 

-623.8616048933 

-683.456632182920 

25 

-609.53231995807574620568 

-677.9637661344 

-743.147260969064 


26 -659.90733322632780520901 -734.3159301916 -805.337898245040 

27 -712.28234551602655145614 -792.9180967274 -870.028542951686 


28 


-766.65735693155709991040 


14 


-853.7702654564 


-937.219194199135 






























TABLE II; Coefficients a'^^\ and ..., 67 ^®^ in the ten-term expansion of the 

E((5, A^e)—function (see, Eq.(12)). 


coefficients 

Ne = 2 

Ne = 3 

Ne = 4: 

Oq 

- 1.0000000000 

0.6249999961 

-0.1576662873 

-1.1249999945 

1.0228047058 

-0.4081459572 

-1.2499999978 

1.5592739061 

-0.8771019520 

bf^'> 

J^e) 

AN.) 

’^3 

,iN.) 

AN.) 

°6 

f^iN.) 

0.0086962674 

-0.0008567981 

-0.0012640360 

0.0003957668 

-0.0030653579 

0.0037323637 

-0.0027229312 

-0.0170061678 

-0.0341594959 

-0.1084597316 

0.2941494433 

-1.4458330710 

2.9184336113 

-3.3508481226 

-0.0429178704 

-0.1694257553 

-0.3362058744 

0.6932000874 

-5.7572871897 

14.3319856750 

-22.3281869042 
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